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We investigate the crossover from Bardeen-Cooper-Schrieffer (BCS) pairing to a Bose-Einstein 
condensate (EEC) in a relativistic superfluid within a boson-fermion model. The model includes, 
besides the fermions, separate bosonic degrees of freedom, accounting for the bosonic nature of 
the Cooper pairs. The crossover is realized by tuning the difference between the boson mass and 
boson chemical potential as a free parameter. The model yields populations of condensed and 
uncondensed bosons as well as gapped and ungapped fermions throughout the crossover region for 
arbitrary temperatures. Moreover, we observe the appearance of antiparticles for sufficiently large 
values of the crossover parameter. As an application, we study pairing of fermions with imbalanced 
^ \ populations. The model can potentially be applied to color superconductivity in dense quark matter 

^— ■>) . at strong couplings. 
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lO . I. INTRODUCTION 

. An arbitrarily weak attractive interaction between fermions in a many-fermion system leads to the formation 
J — ' of Cooper pairs. This phenomenon is well described within Bardeen-Cooper-Schrieffer (BCS) theory [l||. In this 
^\ . situation. Cooper pairs are typically of a size much larger than the mean interparticle distance. The picture changes 
^) • for sufficiently large interaction strengths. In this case. Cooper pairs become bound states, and superfluidity is 
', realized by a Bose-Einstein condensation (BEC) of molecular bosons composed of two fermions. A crossover between 
the weak-coupling BCS regime and the strong-coupling BEC regime is expected 01 • 

Experimentally, this crossover has been studied in systems of cold fermionic atoms in a magnetic trap, where the 
coupHng strength can be tuned around a Feshbach resonance with the help of an external magnetic fleld [3]. Recently, 
these studies have been extended to the case two fermion species with imbalanced populations In this case, the 
crossover is most likely replaced _by one or more phase transitions, and the appearance of exotic superfluids seems to 
Q ' be a very interesting possibility [3, 0] ■ 
^ ' Besides the nonrelativistic atomic systems, there is also a strong motivation to study the relativistic BCS-BEC 
Ch ' crossover. One possible realization is pion condensation, which, for large isospin densities, crosses over into Cooper 
^ I pairing of quarks and antiquarks Another possibility is dense quark matter which may be present in compact stars 
'-j ■ [3| . Under astrophysical conditions of densities of a few times the nuclear ground state density and comparably small 
K> I temperatures of 1 MeV and lower, quark matter is a color superconductor Analogous to electrons in a metal or 

alloy or fermionic atoms in a magnetic trap, quarks form Cooper pairs due to an attractive interaction, here mediated 
by gluon exchange. Because of asymptotic freedom, color superconductivity at asymptotically large densities can 
be studied in a weak-coupHng approach using perturbative methods within QCD [ill, [l3|- However, for moderate 
densities as present in compact stars, the vaHdity of these results is questionable. More phenomenological models such 
as the Nambu-Jona-Lasinio (NJL) model, have therefore been employed, mimicking the gluon exchange by a pointlike 
interaction between the quarks (see Ref. [l3| and references therein). Both QCD and NJL approaches usually are 
applied within a BCS-Hke picture. However, quark matter in compact stars may well be in a strong-coupling regime 
where a BEC-like picture is more appropriate (l3. [isl. \lK \v^\l^. 

In order to describe the crossover from BCS to BEC we shall not consider a purely fermionic model which may 
describe this crossover as a function of the fermionic coupling strength. We rather set up a theory with bosonic and 
fermionic degrees of freedom. Here, fermions and bosons are coupled through a Yul-cawa interaction and required to 
be in chemical equilibrium, 2/i = fib, where fi and fib are the fermion and boson chemical potentials, respectively. 
We treat the (renormalized) boson mass mf, and the boson-fermion coupling g as free parameters. Then, tuning 
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the parameter x = —{ml^ ~ 1^1)1(^9^) drives the system from the BCS to the BEC regime. The fermionic chemical 
potential shall be self-consistently determined from the gap equation and charge conservation. This picture is inspired 
by the boson-fermion model of superconductivity considered in Ref. [l^ , which has been used in the context of cold 
fermionic atoms [23| • It also has possible applications for high-temperature superconductivity [2l| . For simplicity, we 
shall restrict ourselves to the evaluation of the model in a mean-field approximation. 

The paper is organized as follows. In Sec. HIl we set up the simplest version of the model, taking into account a 
single fermion species. We define the crossover parameter in Sec. Ill CI and derive the density and gap equations in 
Sec. Ill Dl The solutions of these equations are presented in Sec. IIIII We consider a vanishing temperature in Sec. 
IIII A[ present the crossover at the critical temperature in Sec. IIIIB^ and show results for a fixed crossover parameter 
and arbitrary temperature in Sec. IIII Cl In Sec. IIII Dl we present the results for the ratios A/^ and Tc/A. Finally, 
we extend our model to two fermion species in Sec. HVl This extension allows us to consider pairing of fermions with 
imbalanced populations, which is an unavoidable complication in quark matter at moderate densities (23 |. 

Our convention for the metric tensor is g^'^ = diag(l, — 1, — 1, — 1). Our units are S = c = fcs = 1. Four- 
vectors are denoted by capital letters, K = K^^ = (fco,k) with k = |k|. Fermionic Matsubara frequencies are 
ujn = iko = {2n + 1)ttT, while bosonic ones are a;„ = ika = 2mrT with the temperature T and n an integer. 



II. THE BOSON-FERMION MODEL FOR A RELATIVISTIC SUPERFLUID 

A. Setting up the model 

We use a model of fermions and composite bosons coupled to each other by a Yukawa interaction. The Lagrangian 
is given by a free fermion part Cf, a free boson part Cb and an interaction part £/, 

C = Cf + Cb + Ci, (1) 

with 

£f = ilj{i'-f^df, + 7oAf - m)ip , (2a) 

Cb = \{dt-iHb)^\^ -\\7^\^ -mtl^W (2b) 

Ci = g{ipipcij5-ip + ip*i>ij5ipc) ■ (2c) 



The fermions are described by the spinor ip, while the bosons are given by the complex scalar field ip. The charge 

'j-< 

conjugate spinors are defined by tpc = dtp and ipc ~ V'^C' with C = 17^7°. The fermion (boson) mass is denoted 
by m (mb). We choose the boson chemical potential to be twice the fermion chemical potential. 



fib = 2fi. (3) 

Therefore, the system is in chemical equilibrium with respect to the conversion of two fermions into one boson and 
vice versa. This allows us to model the transition from weakly-coupled Cooper pairs made of two fermions into a 
molecular difermionic bound state, described as a boson. The interaction term accounts for a local interaction between 
fermions and bosons with coupling constant g. In order to describe BEC of the bosons, we have to separate the zero 
mode of the field if Moreover, we shall replace this zero-mode by its expectation value 

cP = (^0) (4) 

and neglect the interaction between the fermions and the non-zero boson modes. This corresponds to the mean-field 
approximation. Then, with the Nambu-Gorkov spinors 

*=f^t')' * = (^'^c), (5) 



'c , 

the Lagrangian can be written as 

C = ^-^S-'^ + [^,1 - mlM' + \idt- i^ib)v\^ - |V^|2 - ml\ipf . (6) 

Note that we have dropped the mixing terms of zero and non-zero boson modes since they vanish when carrying out 
the path intergal. Here is the inverse fermion propagator which reads in momentum space 
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It is instructive to compare this form of the propagator to the corresponding one in a purely fermionic model, see for 
instance Ref. As expected, the Bose condensate is related to the diquark condensate A, 



A = 25(/) . 



(8) 



As we shall see below, cf. Eq. I|lip . A is the energy gap in the quasi-fermion excitation spectrum. In a purely fermionic 
model, A = 2G{ipcil5''P)i where G is the coupling constant related to the interaction between the fermions. Note 
that G has mass dimension —2, while our boson-fermion coupHng g is dimensionless. Therefore, g does not play the 
role of the crossover parameter, as G does in the fermionic model. We shall explain this in more detail in Sec. Ill CI 



B. Thermodynamic potential 



In order to obtain the thermodynamical potential density fl, we compute the partition function 

»1/T 



[(m'][d^][dip][dip*]cxp 



drd-^xC 



(9) 



where T is the temperature, and C is the Lagrangian in the mean field approximation given in Eq. ([6]). The thermo- 
dynamic potential density is then obtained from Q. = —T/V InZ, where V is the volume of the system. One obtains 
after performing the path integral and the sum over Matsubara frequencies. 



^ = -E 



e=± 



d-^k 

(2^ 



2Tln 



l+cxp( 



(ml -t^DA^ 
4g2 



-T 



d^k 
(2^ 



2rin 



1-exp ( 



(10) 



We have used Eq. |[8]) and denoted the quasi-particle energy for fermions (e — +1) and antifermions (e = — 1) by 



and the (anti) boson energy by 



= + ml - efit ■ 
Furthermore, we have assumed A (and thus 0) to be real. 



(11) 



(12) 



C. Crossover parameter 

We shall now define the crossover parameter whose variation carries the system from the BCS to the EEC regime. 
Let us first recall the corresponding crossover parameter in a purely fermionic model. In this case, the fermionic 
coupling G has to be renormalized. This is in contrast to the weak-coupling regime where the gap equation is well- 
defined with the bare coupling G (a natural cutoff is provided by the Debye frequency in the non-relativistic case; in 
QCD, the gap is a function of momentum and peaks around the Fermi surface, providing a regular behavior of the 
gap equation). The renormaHzed coupHng is proportional to the scattering length. In the context of cold fermionic 
atoms, the scattering length is the physical quantity which can be controlled upon tuning the external magnetic field. 
For the relativistic case, see Ref. [l7| for the relation between the renormalized coupling and the scattering length. 
The definition of the crossover parameter in the present model goes along the same lines. Instead of a renormalized 
coupling we introduce the renormalized boson mass 



'fc.r 



4.9 
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A=/j=T=0 

This allows us to define the (renormaHzed) crossover parameter 



(13) 



(14) 



4 



The parameter x can be varied from negative values with large modulus (BCS) to large positive values (BEC). In 
between, a; = is the unitary limit [13, HI, H^- Therefore, 1/x behaves similar to the scattering length: the BCS 
(BEC) Hmit is approached via 1/x t {1/x [ 0) while the unitary regime corresponds to 1/x = ±oo. 

We may thus write the thermodynamical potential in terms of the parameters {x,g) instead of the original pair 
(mh, g). To this end, we have to express the bare boson mass rrib in terms of x and g. With the help of Eqs. ifTSj) and 
(fT4)) we find 

ml - Mb = 4.92 [j -x \ = Ag^xo - x) . (15) 

For sufficiently small fermion masses, to ^ A, where A is the cutoff in the momentum integrals, we have 

xo ^ ^ . (16) 

One sees from Eq. ifTSj) that xq is an upper Hmit for x in order to ensure non-negative bosonic occupation numbers. 
Moreover, in the limit of large x — > xg the boson chemical potential approaches the (bare) boson mass, — > "T-b- The 
condition of Bose-Einstein condensation in a free bosonic system with fixed bosonic charge is fib = "mb- In the present 
model, however, we shall observe a nonzero Bose condensate also for fib < mb, corresponding to x < xq- 

Having defined the crossover parameter x and its definition range x S [— oo,a;o], we note that, within our simple 
model, we are left with the second free parameter g. We shall discuss below how the choice of g effects the behavior 
of the system in the BCS-BEC crossover. For most of our results we shall, however, use a single fixed value of g. 



D. Densities and gap equation 



Next, we derive the charge conservation equation and the gap equation which shall later be solved numerically. The 
total charge density 

— f 

can, using Eq. fTO|) . be written as 

n = np + no + ub ■ (18) 

Here, the fermionic contribution is given by 

(^^[/f(4)-/f(-4)], (19) 

where we abbreviated 

C^ = efco-e/i, (20) 

and is the Fermi distribution function, /i?(a;) = l/[cxp(x/r) + 1]. The factor 2 in front of the sum in Eq. lfT9|) 
originates from the two spin degrees of freedom. From Eq. I|19p one recovers the limit case of a free Fermi gas at zero 
temperature, 

nF(A = T = 0) = ^ [e(M - to) - e(-M - "^)] • (21) 



The condensate density is 



and the thermal boson contribution is 



"0 = — 5— , (22) 
5 



nB^2Y^e I j^IbH), (23) 
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where /b is the Bose distribution function, fsix) = l/[exp{x/T) — 1]. The factor 2 in the boson densities originates 
from Eq. ([3]), i.e., from the fact that each boson is composed of two fermions. For the following, let us also define the 
charge fractions 

- "-s/F _ no .^ 

Pb/f = , po = —, 24 

' n n 



and an effective Fermi momentum pp through 



(25) 



The effective Fermi energy is then given by = \/Pf + The various densities appearing on the right-hand 
side of Eq. lfT8|) are interpreted as follows. The fermions that contribute to np are, for temperatures below the 
superfiuid transition temperature, constituents of weakly coupled Cooper pairs. For temperatures larger than the 
transition temperature, np corresponds to free fermions. The bosons that contribute to the boson density Hb are, for 
all temperatures, uncondensed molecular bound states, composed of two fermions. Condensation of these pairs can 
only occur below the transition temperature and results in a nonzero condensate density uq. 

In order to find the gap equation in the case of a fixed charge density, we have to minimize the free energy density 

F = n + ^m. (26) 

Here, is an implicit function of n (and of the gap A) through Eq. lfT7|) . Minimization with respect to A yields 

dA ~ 9A ^ 9^aA ^""M ~ aA ' ^ 

where Eq. |[T7| has been used. For A 7^ Eq. ([271 reads 



E 

e=± 



TT^tanh-A,-— . (28) 



(2^)3 \2el 2T 2eko 



Note that the density n in Eq. ifTSj) was obtained by taking the derivative with respect to /i at fixed nih (not at fixed 
x). This is necessary to obtain a nonzero condensate contribution uq. Also the equivalence of dF/dA and dfl/dA in 
Eq. l(27|) is obtained under this premise. At fixed nib we get /x by solving Eqs. (fT8| and l|28p . We can then obtain x 
from mi, and fj. via Eq. lfT4|) . In this way we have a one-to-one mapping between rri}, and x. Here we should emphasize 
that the current case is different from that x is fixed from the very beginning before Eq. (fT8|l is derived. 



III. RESULTS AND DISCUSSION 



The two coupled equations (flSl) and l(28| with the definitions (fTOl) . l(22|l . and i(23l) shall be used in the following to 
determine the gap A, and the chemical potential fi as functions of the crossover parameter x, see Eq. (fT4| . and the 
temperature T at fixed effective Fermi momentum pp, fermion mass to, and boson-fermion coupling g. The solution 
[A(x,T), fi{x,T)] can then, in turn, be used to compute the densities of fermions and bosons in the x-T plane. We 
shall present results for the zero-temperature case, Sec. lIII Al and at the critical temperature Tc, Sec. IIII Bl Then, we 
show results for a fixed x and arbitrary temperature T, Sec. IIII CI Throughout these subsections, we shall fix 

^=0.3, "^=0.2, ff = 4. (29) 

In Sec. IIII P] we present the ratios Aq/^q and Tc/ Aq for different values of m and g as a function of x. Here and in 
the following, we use the subscript at A and yU to denote the zero-temperature value. 

A. Zero temperature 

For T = 0, there are no thermal bosons, ub — 0, and Eqs. ifTSj) and ([28l) become 

n = np + riQ , (30a) 




(30b) 
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with the zero-temperature expressions for the fermion densities 



np — 




(31) 



The numerical results for the solution of the coupled equations pOa|) and i|30bp are shown in Fig. [TJ The left panel 
shows the fermion chemical potential ^0 and the gap Aq as functions of a;. In the weak-coupHng regime (small x) we 
see that the chemical potential is given by the Fermi energy, /io = e_F- For the given parameters, cf/A ~ 0.36. The 
chemical potential decreases with increasing x and approaches zero in the far BEC region. The gap is exponentially 
small in the weak-coupling region, as expected from BCS theory. It becomes of the order of the chemical potential 
around the unitary limit, x = 0, and further increases monotonically for positive x. In the unitary limit, we have 
li/ep — 0.37, while in nonrelativistic fermionic models /i/e_F — 0.4 — 0.5 was obtained [2i. [25l. [2^ . 

The corresponding fermion and boson densities are shown in the right panel of Fig. [TJ These two curves show 
the crossover: at small x all Cooper pairs are resonant states, which is characterized by a purely fermionic density, 
n = np; at large x, on the other hand. Cooper pairs are bound states and hence there are no fermions in the system. 
The charge density is rather dominated by a bosonic condensate, n = uq. The crossover region is located around 
X = 0. We can characterize this region quantitatively as follows. We write the boson mass as mi, = 2m — i?bind- 
Then, a bound state appears for positive values of the binding energy i?bind, i.e., for 2m > mb. With Eq. (flSl) this 
inequality reads 



Since /iq is a monotonically decreasing function of x, this relation suggests (note that xq — x > by construction): 
(i) for sufficiently large x < xo bound states appear for any fixed g {ii) the larger g the "later" (= larger values of x) 
bosonic states appear. We have confirmed these two statements numerically by using different values of g. The value 
(7 = 4 is chosen such that there is an approximately balanced coexistence of fermions and bosons aX x ~ Q, np ~ no, 
as can be seen in the right panel of Fig. [TJ 

One may ask whether there is a contribution of antifermions to the total fermion charge. In the BCS regime 
there is a Fermi surface given by ^ > and antifermion excitations are obviously suppressed. However, during the 
crossover, ^ decreases and there might be the possibility of the appearance of antifermions. The contributions of 
fermions and antifermions to the total fermion charge seem to be given by the terms e = +1 and e = — 1 in Eq. (|3T|) . 
A separate discussion of these terms is not straightforward because they contain divergent contributions which cancel 
in the sum but not in each term separately. Thus, a renormalization of both terms would be required. In the BCS 
regime, a; —00, vacuum contributions cx A'^ have to be subtracted. For nonzero values of the gap, however, more 
divergent terms appear involving powers of both the cutoff and the gap. (Note that this problem is not unlike the one 
encountered in Ref. [23| , where medium-dependent counter terms were introduced in the calculation of the Meissner 
mass. In fact, we shall choose a similar renormalization in the following calculation of the energy density.) 

In any case, separate charges of fermions and antifermions are not measurable quantities since any potential experi- 
ment would solely measure the total charge. Therefore, we shall describe the onset of a nonzero antiparticle population 
in terms of the energy density. In this quantity, we expect the contributions from particles and antiparticles to add 
up, in contrast to the charge density where the contributions (partially) cancel each other. 

Let us first discuss the quasi-fermion and quasi-antifermion excitation energies given by and from Eq. (fTTj) . 
Inserting the numerical solutions for ^ and A into these energies results in the curves shown in Fig. [2j These 
excitation energies show that, for large values of x, quasi-fermions and quasi-antifermions become degenerate due 
to the vanishingly small chemical potential. Because of the large energy gap, we expect neither quasi-fermions nor 
quasi-antifermions to be present in the system. 

This statement can be made more precise and generalized to nonzero temperatures upon considering the energy 
density E. Using the thermodynamic potential density 17 from Eq. (fTOj) and the entropy density S = —dn/{dT) we 
have E = fl + /in + TS. We obtain 



^il < m^ -g^{xo - x). 



(32) 



E — Ep + Eb , 



(33) 



with the fermionic and bosonic contributions 




(34a) 



Eb 




(2^ 



c.^,[l + 2/B(t^D] + (^o-x)A2 



+ /i(no + ub) ■ 



(34b) 
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The renormalization of the fermionic part can now be chosen such that there are no quasi-particles at T = np = 0, in 
accordance with the above argument. Hence we subtract the "vacuum contribution" Ep{T = np = 0) to obtain the 
renormaUzed energy density 



Ef ». — 2 



(27r)3 



lfF{4) + ^J-np. 



(35) 



For T = 0, only the second term survives (remember that > 0), and Ep behaves as shown in the left panel of Fig. 
m For nonzero temperatures, however, we see that there is a nonzero fermionic energy density even for np = 0. This 
is related to the excitation of quasi-antifermions, as we shall discuss in the next subsection. 

For the bosonic energy density, we subtract the analogous vacuum part Eb{T — no + ns — 0). Hence we obtain 
the renormalized energy density 



Ef 



(2^ 



^fc/B(^fe)+M('^o + ni3). 



(36) 



At r = 0, we have Eb.t = ^'^■0, which is shown in the left panel of Fig. [H We see that the bosonic energy density 
vanishes in the BCS regime because there is no Bose condensate in this case, n^ = 0. In the far EEC regime, where 
no 7^ 0, the energy density vanishes too because the boson chemical potential, since coupled to the fermion chemical 
potential, vanishes. Only in the crossover region, where both the condensate and the chemical potential are nonzero, 
the energy density is nonvanishing. 





Figure 1: (Color online) Crossover at zero temperature from the BCS regime (small x) to the BEC regime (large x). Left 
panel: fermion chemical potential /io (blue dotted) and gap Ao (red dashed) in units of effective Fermi energy ep- Right panel: 
condensate fraction (red solid), fermion fraction (blue solid). 



B. Critical temperature 

In this section, we calculate the critical temperature Tc and the corresponding particle densities as functions of x. 
Upon setting A = in the charge density equation IjlSp and the gap equation l(28|) one obtains 



e=± 

with the fermion density 



n = np + iiB , (37a) 
(2^)3 \2^l 2T, 2eko, 



V /T^f^tanhiL-^), (37b) 



"^-25:e /^/p(CD, (38) 
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Figure 2: (Color online) Fermion and antifermion excitation energies and as defined in Eq. I|lip for three different values 
of the crossover parameter x/xo at T = as a function of the momentum k (both and k are given in units of ef). In the BCS 
regime (left panel) the energy gap is small and the fermion excitations are well separated from antifermion excitations. Both 
excitations approach each other in the unitary regime (middle panel), and become indistinguishable in the far EEC regime 
(right panel). Note in particular that the minimum of the antiparticle excitation is not a monotonic function of x. 
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Figure 3: (Color online) Energy density of fermions and bosons in units of epn at T = (left panel) and at T = Tc (right 
panel). In the BCS regime, Ep ~ epn and Eb = for all temperatures. The large fermionic and bosonic energy densities in 
the EEC regime at T = Tc indicate the occupation of (quasi-)antiparticle modes. 



and the boson density given by Eq. l(23|) . Strictly speaking, the original gap equation l(28|) is only valid for nonzero 
A (in its derivation, one has to divide by A). Therefore, Eq. (|37bp has to be understood as a limit for approaching 
the critical temperature from below, T "f Tc, i.e., for infinitesimally small A. Eqs. I|37ap and ()37bp can now be used 
to determine Tc and the corresponding chemical potential ^c- 

The results are shown in the left panel of Fig. HI We see that the chemical potential behaves qualitatively as for 
zero temperature. The critical temperature, while exponentially small in the BCS regime, becomes of the order of 
and then larger than the chemical potential during the crossover. This is one of the characteristics of the strong 
coupling regime and one reason why this model (in its nonrelativistic version) is used to describe high-temperature 
superconductivity [IH- In Sec. lIIIDl we use the ratio Tc/ Aq to illustrate the high-Tc behavior. 

The right panel of Fig. |4] shows the particle density fractions for fermions and bosons. A crossover similar to the 
zero-temperature case can be seen. The density fractions of fermions and bosons suggest that the crossover is shifted 
to a slightly larger value of x compared to the zero-temperature case. While at zero temperature np = m occurs at 
x/xo ~ 0, here we have np = ns at x/xo — 0.3. It is clear that there is no Bose condensate at T = Tc; the bosonic 
population rather consists of thermal molecules. These are uncondensed, strongly-coupled Cooper pairs (see Ref. 
for a recent discussion of this effect in the context of cold atoms) . We see that uncondensed pairs do not exist in the 
BCS limit. In this case, the superfluid phase transition occurs "abruptly", with pair formation and condensation at 
the same temperature. 

In the right panel of Fig. [3] we show the fermion and boson energy densities. These curves are obtained by inserting 
the solutions for fic and Tc into Eq. (|35|) and l(36|) and making use of A = 0. We see that, in contrast to the zero- 
temperature case, the fermionic energy density increases with x despite np —* 0. This is easy to understand from the 
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corresponding quasi-fermion excitations. For all temperatures, the vanishing chemical potential renders quasi-particles 
and quasi-antiparticles degenerate. Whereas at T = they are both gapped by A, at T = Tc they are gapped only 
by the fermion mass m. For large values of x we have Tc :$> m and quasi-fermions as well as quasi-antifermions can 
be thermally excited. In this context, it would be interesting to consider the formation of chiral condensates which 
may be initiated by the degeneracy of quasi-fermions and quasi-antifermions. We leave this extension of the model 
for future studies. 

The large increase of the bosonic energy density can be understood in the same way. Note, however, that, in contrast 
to the fermion mass, the boson mass decreases with increasing crossover parameter x. This difference, together with 
the different statistics of bosons and fermions gives rise to the qualitatively different behavior of Eb.t compared to 
Ep r. The strong increase of antiparticle densities has also been predicted in other models and has been termed 
"relativistic BEC (RBEC)" The relativistic effects have also been studied in Ref. 





-0.2 



Figure 4: (Color online) Crossover at the critical temperature. Left panel: fermion chemical potential /ic (blue dotted) and 
critical temperature Tc (red dashed) in units of the effective Fermi energy ef. Right panel: fermion fraction (blue solid), 
thermal boson fraction (red solid). 



C. Fixed coupling 

In the previous two subsections we have presented the solution of Eqs. flS]) and (|28|) along two lines in the x-T 
plane: along the line T = (Sec. lIII A|) and along the (curved) line T = Tc (Sec. IIII B|) . Now we explore a third path 
by fixing the crossover parameter and vary the temperature from zero to values beyond Tc. We shall use x/xq = 0.2 
which is in the intermediate-coupling regime, where both fermionic and bosonic populations are present. For T < Tc, 
we use Eqs. ifTS)) and l[28|) to determine fi and A. For T > Tc, the gap vanishes, A = 0, i.e., we have the single 
equation l|37ap to determine the chemical potential fi. 

The condensate and the fermion and boson density fractions are shown in Fig. [H At the left end, T = 0, one 
recovers the results shown in Fig. [T] at the particular value x/xg = 0.2, while the point T/Tc = 1 reproduces the 
respective result shown in Fig. IH The second-order phase transition manifests itself in a kink in the density fractions 
and a vanishing condensate. Below Tc we observe coexistence of condensed bound states, condensed resonant states, 
and, for sufficiently large temperatures, uncondensed bound states. We obtain thermal bosons even above the phase 
transition. They can be interpreted as "preformed" pairs, just as the uncondensed pairs below Tc. This phenomenon 
is also called "pseudogap" in the literature [H, [2§|- It suggests that there is a temperature T*{x) which marks the 
onset of pair formation. This temperature is not necessarily identical to Tc. In the BCS regime, T*{x) = Tc{x), while 
for a; > 0, T*(x) > Tc{x). Of course, our model does not predict any quantitative value for T* because thermal bosons 
are present for all temperatures. Therefore, we expect the model to be valid only for a Hmited temperature range 
above Tc. 



D. The ratios Ao//io and Tc/ Ao 



We finally present the results for the ratios Aq/^xq, and Tc/ Aq. They shall serve as a discussion of the dependence 
of our results on the boson-fermion coupling g and the fermion mass m. Both g and m were fixed throughout the 
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Figure 5: (Color online) Density fractions in the crossover regime at fixed xjxa — 0.2 as functions of temperature: condensed 
bosons (red solid), fermions and uncondensed bosons (blue dotted and red dashed, respectively). 



previous sections. Moreover, we shall see that we reproduce values of these ratios obtained in different models in 
certain limit cases. 

Fig. [6] shows the ratio Ao//io, using the results for Ag and /io from Sec. lIII Al From both panels one can read off the 
value of the ratio in the unitary limit, a; ^ 0. For the fermion mass that has been used in the previous subsections, 
to/A ~ 0.2, we find 1.2 < Ao//io ^ 1-4. The exact value depends on the choice of g. This range is in agreement 
with nonrelativistic, purely fermionic models [13, 113, HBl • The right panel shows that the ratio in the unitary limit 
decreases with decreasing fermion mass. In particular, in the ultrarelativistic limit to = we find Aq/^q — 1-0. 

In Fig. [7] we show the ratio Tc/Aq, using the result for Tc from Sec. IIIIBl From BCS theory we know 

lim ^ = — ~ 0.57 , (39) 

a;^-oo Aq TT 

where 7 ~ 0.577 is the Euler-Mascheroni constant. This value is reproduced in our results, independent of g and to. 
Upon increasing the crossover parameter x, the ratio deviates from its BCS value and increases substantially during 
the crossover where it strongly depends on the coupling g. Therefore we make no predictions for its value in the 
unitary regime. However, we see that in the EEC regime, the value again becomes independent of the parameters 
and assumes a value 

lim ^ ~ 0.50 . (40) 

x^xq Ag 



IV. TWO-SPECIES SYSTEM WITH IMBALANCED POPULATIONS 



It is straightforward to extend our boson-fermion model to two fermion species with cross-species pairing. This 
allows us to introduce a mismatch in fermion numbers and chemical potentials which imposes a stress on the pairing. 
This kind of stressed pairing takes place in a variety of real systems. For example, quark matter in a compact star is 
unlikely to exhibit standard BCS pairing in the color-fiavor locked (CFL) phase, i.e., pairing of quarks at a common 
Fermi surface. The cross-fiavor (and cross-color) pairing pattern of the CFL phase rather suffers a mismatch in 
chemical potentials in the pairing sectors hu — rs and bd — gs (r, g, b meaning red, green, blue, and u, d, s meaning 
up, down, strange). This mismatch is induced by the explicit fiavor symmetry breaking through the heaviness of 
the strange quark and by the conditions of color and electric neutrality. Our system shall only be an idealized and 
simplified model of this compHcated scenario. However, as in the previous sections, we shall allow for arbitrary values 
of the crossover parameter and thus model the strong coupling regime of quark matter. We shall fix the overall charge 
and the difference in the two charges. This is comparable to the effect of neutrality conditions for matter inside a 
compact star, which also impose constraints on the various color and flavor densities. Our focus will be to flnd stable 
homogeneous superfluids in the crossover region and, by discarding the unstable solutions, identify parameter values 
where the crossover in fact becomes a phase transition. We shall restrict ourselves to zero temperature and defer the 
full analysis of the two-species system to a future study [sOl • 
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Figure 6: (Color online) Ratio of gap over chemical potential at zero temperature for crossover parameters 2;/xo < on a 
logarithmic scale. The left end of the horizontal axis corresponds to the BCS regime, the right end, where a; — » 0, corresponds 
to the unitary regime. Left panel: ratio for different values of the boson-fermion coupling g. Right panel: ratio for different 
values of the fermion mass. 
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Figure 7: (Color online) Ratio of critical temperature over zero-temperature gap throughout the BCS-BEC crossover. Left 
panel: ratio for different values of the boson-fermion coupling g. Right panel: ratio for different values of the fermion mass. 



A. Two-fermion system 
We start by replacing the fermion spinor ip in the Lagrangian ^ with a two-component spinor 



and the chemical potential /i with the matrix diag(/ii, ^2)- The factor 1/V2 accounts for the same normalization for 
the total fermion number as in the case of single fermion species. We assume both fermions to have the same mass m. 
Cross-species pairing is taken into account in the interaction part of the Lagrangian £/ in ([2c|) which is now replaced 
by 

Ci = giiptpchbcri'ip + V'*tp-i75cri'4'c) , (42) 

where the Pauli matrix ai is a matrix in the two-species space. We denote the average chemical potential and the 
mismatch in chemical potentials by 

M = , Sfi = . (43) 
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Then, the bosonic chemical potential is 

Mb = 2/1 . 

The thermodynamic potential differs from the one-fermion case in the dispersion relation for the fermions, 



(44) 



n 



4,92 ^ 2 



wf, + 2rin 1-e 



(27r)3 

1 4. e-(^l+^^)/^ 



Tin 



1 + e" 



-(e|-5,x)/T 



}• 



where 



and Lo^ as in Eq. I|12p . At zero temperature f2 becomes 



4^2 



(27r)3 



[e^fc + (<5/i-e^)e(<5/i-4)] 



The particle number densities for each species are derived from the thermodynamic potential, 

dn no 



where no is given by Eq. (|22|) and 

"/,l/2 



(2^ 



2el 



2ef 



We shall evaluate the model for fixed sum and difference of the particle number densities 

d^k 



n = rii + 712 = - 
5n = ni — 712 = — 



dSfi 



3 ie(4-5/.). 



d^k 



e(d> - el) . 



(45) 



(46) 



(47) 



(49) 

(50a) 
(50b) 



Of course, the bosons contribute equally to both particle numbers and thus do not appear in Sn. The gap equation 
becomes 



E 



d^k 
(2^ 



1 1 



2ef. 2efco 



(51) 



We shall solve Eqs. ((50| and fSTj) for the variables /j,, (5/x, and A. 



B. Possible Fermi surface topologies 

Before we solve the equations, let us comment on their structure, in particular the appearance of the step function. 
It is convenient to rewrite the step functions such that their effect can be translated into the boundaries of the dk 
integration. Furthermore, it is instructive to discuss the different particle and antiparticle occupation numbers in 
momentum space with the help of the step functions. As we shall see, the occupation numbers are discontinuous at 
the zeros of the dispersion relation — 5^. We assume without loss of generality that J1,S^ > (in the numerical 
solution we ensure this by choosing n,Sn > 0). We abbreviate 



(52) 
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cases 


characteristics 


parameter region 


NJ^i{k,) (lermionsj 


Nf -{k) (anti-fermions) 


I 


fully gapped 


Sjj. < A 01 C,+ < m 


iV+ = iV+ = N+, 


iV-, = Ny, = iVg-p 


II 


breached for fermions 
gapped for anti-fermions 


m < C- < C+ 


N+ = 1, iv+ = 0, G [ypT, ypT] 
iv+ = N% = iv+p, k i [ypT, yp^] 


same as case I 


III 


single EFS for fermions 
gapped for anti-fermions 


IC-I < m < C+ 


iV+ = 1, iV+ = 0, fc G [0, VpT] 


same as case I 


IV 


single EFS for fermions 
single EFS for anti-fermions 


m < -C- < C+ 


same as case III 


iv^, = 0, iV72 = -i, fcG[o,ypr] 

= = ^^gap, k i [0, ypT] 



Table I: Four possible parameter configurations and corresponding fermion (e = -I-) and antifermion (e — — ) occupation 
numbers Nf i{k) {i — 1,2). The occupation numbers are defined as the integrand in Eq. l|49|l . We abbreviate "effective Fermi 
surface" by EFS and N^^p = e{e% — £_l)/{2e%). The different cases are illustrated in Fig. [8] 



Then, the step functions are 

9(5/. -6+) = e(5M-A)[e(p+)e(VpT-fc)-e(C-)e(p_)e(VP^-fc)] , (53a) 

Q{5^l~eZ) = Q{5^i~/\)Q{-Q-)Q{p-)Q{^-k). (53b) 

We see that the step functions only give a contribution if > A. Three different terms occur, each corresponding to a 
different scenario, distinguished by the topology of the effective Fermi surfaces of fermions and antifermions. Together 
with the fully gapped state, these are four possible cases. We list these cases and their characteristics in Table [J and 
Fig. [HI The fermion dispersion (e = -I-) has either zero, one, or two zeros. A zero corresponds to an effective Fermi 
sphere. In particular, a fully gapped state is characterized by the disappearance of an y Fe rmi surface. The case of two 
effective Fermi surfaces is termed "breached pairing", following the usual terminology [31|. The antifermion dispersion 
(e = — ), in contrast, can have either zero or one effective Fermi surfaces. Here, the asymmetry between fermions and 
antifermions is given by the choice Jl > 0. Hence there is no breached pairing for antifermions. We find the interesting 
possibility of filled Fermi surfaces for fermions of species 1 and antifermions of species 2, see case IV in Tableland 
lower right panel in Fig. [H In Sec. IIV Dl we shall see that this is indeed a stable solution for certain values of the 
crossover parameter. 



C. Number susceptibilities 

In order to check the gapless states for their stability, we have to compute the number susceptibility matrix [1, [3| 

X.-^, . = 1,2. (54) 

Note that we fix n and Sn (or equivalently ni and n2) in our solution. Hence x can be regarded as measuring the 
response of the system to a small perturbation away from this solution. In particular, a stable solution requires the 
mismatch in density to increase for an increasing mismatch in chemical potentials. Therefore, a negative eigenvalue 
of this 2x2 matrix indicates the instability of a given solution. The susceptibility is evaluated upon using 

- lc?no duf^i _ / 2/lA dufA dA duf^j 

~ 2 d^lJ dnj ~^ 2g^ \ dA J dfij dfij ' ^ ' 

The partial derivatives of the fermion densities with respect to the gap and chemical potentials are straightforwardly 
computed with the help of Eq. l(49|) . The partial derivative of the gap with respect to the chemical potentials can be 
computed from the gap equation. To this end, we take the (total) derivative with respect to /ii/2 on both sides of Eq. 
([STI and solve the equation for dA/d^i/2- The results for the various terms are given in Appendix [Al 



D. Stable and unstable gapless superfluids 



We first present the results for a certain mismatch in densities. 
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Figure 8: (Color online) Occupation numbers for fermions (red) and anti-fermions (blue). Solid and dotted lines are for species 
1 and 2 respectively. The panels illustrate the qualitatively different cases I through IV from Table [H Note that dotted and 
solid lines coincide in various regions, e.g., for all k in the upper left panel. See last two columns of TableUfor the precise form 
of the occupation numbers. 



Moreover, we use n = p^/{3Tr'^), and pp, m, g, asin the main part of the paper, see Eq. l(29|). Then, for T ~ we have 
a set of coupled equations l[50|) and l(5T|) for any value of the crossover parameter x. We solve this set of equations 
for Ji, 6fi, and A. The results are shown in Fig. [9l The left panel shows that the behavior of the average chemical 
potential and the gap are not unlike the case with a single fermion species, cf. Fig. [TJ Moreover, we see that 6^j. > A 
for all X. This is clear from Eq. I|50bp : any nonzero Sn goes along with S^j. > A. In other words, the standard fully 
gapped pairing does not allow (at T = 0) for a difference in fermion numbers. Thus all solutions correspond to gapless 
pairing and case I in Table [J and Fig. [8] does not appear. We have indicated in Fig. [9] for which values of x which of 
the cases II, III, and IV occurs. We have also marked the onset of instability. Evaluation of the susceptibility matrix 
xi'pjS^, A) yields negative eigenvalues for x < —0.04x0 = X- and x > 0.28 xq = x^. In fact, one of the eigenvalue 
diverges at these points. Denoting the two eigenvalues of X by xi, X2, we have 

{+00 for X I x^ , a; I a;+ 
—00 for X 1 X- , X I X-^ 

X2 > for all X . 

Unstable regions of negative xi are shaded in both panels of the figure. We see that the breached pair solution is 
always unstable. This is expected from similar results from mean-field studies for nonrelativistic systems [H] as well as 
quark matter (gapless CFL and 2SC phases [121) ■ Interestingly, the stable region consists of two qualitatively distinct 
states, labelled by III and IV. The system accounts for a given difference in number densities not only by filling an 
effective Fermi surface with particles of species 1 (state III) , but also by additionally filling an anti- Fermi surface of 
species 2 (state IV). An interesting feature of this state is that in the limit of equal Fermi surfaces (1 and anti-2) 
all charges (ni and 712) are confined in the Fermi sphere, just as in the unpaired phase. This can be seen from the 
lower right panel in Fig. [8l the total occupation numbers for momenta larger than the effective Fermi momentum 
vanish since particle and antiparticle contributions cancel each other. Before this limit is reached, however, this state 
becomes unstable for x > x+. This instability in the BEC regime is in contrast to nonrelativistic systems, where a 
gapless solution (there with a single effective Fermi surface) persists throughout the BEC region 



(57a) 
(57b) 
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Figure 9: (Color online) The crossover with imbalanced population at zero temperature. Left panel: JI (blue dashed), Sjj. (black 
solid) and A (red dotted) in units of ep for fixed n and 6n. Right panel: bosonic (red dotted) and fermionic (blue dashed) 
density fractions. The shaded parameter regions are unstable with respect to a negative susceptibility. The states II, III, and 
IV correspond to the respective Fermi surface topologies discussed in Table HI 



The right panel of Fig. [9] shows the bosonic and fermionic number density fractions, cf. Fig. [T] for the analogous 
curves without mismatch. Taking the value of x where po = pp as an indicator, we see that the BCS-BEC crossover is 
shifted to a larger value of a; (x ~ vs. a; ~ 0.1 xq). Of course, what was a crossover in Fig. [1] is now actually replaced 
by at least two phase transitions. The regions marked as unstable will be replaced by a different phase. It is beyond 
the scope of this paper to determine these phases, but it can be expected that they are spatially inhomogeneous, 
for instance a mixed phase where a superfluid and normal phases are spatially separated, or some kind of "Larkin- 
Ovchinnikov-Fulde-Ferrell" (LOFF) state [s^. In the stable region we see that the change in Fermi surface topologies 
(from state III to IV) is visible in the boson and fermion fractions which exhibit a kink at this point. 

We finally present a phase diagram in Fig. [lO] for arbitrary (positive) values of 5n/n. Since we do not consider 
spatially inhomogeneous phases, this phase diagram is incomplete. Its main point is to identify regions where homo- 
geneous gapless superfluids may exist. We find that for sufficiently large mismatches, 5n/n > 0.02 there is a region 
where no solution with nonzero A can be found (we have not shown this region in the above results for 5n/n = 0.5). 
We see that the region of stable superfluids shrinks with increasing Sn/n, as expected. Note that the horizontal axis 
Sn/n = is not continuously connected to the rest of the phase diagram. For vanishing mismatch in densities a stable, 
fully gapped superfluid exists for all a;, as we saw in the main part of the paper. One should thus not be misled by 
the instability for arbitrarily small mismatches. 

We conclude with emphasizing the two main qualitative differences to analogous phase diagrams in nonrelativistic 
systems: (i) within the stable region of homogeneous gapless superfluids there is a curve that separates two different 
Fermi surface topologies; this is the right dashed-dotted line in Fig.[TOl (m) for large x the gapless superfluid becomes 
unstable even in the far BEC region; this is the shaded area on the right side in Fig. [lOl 



V. SUMMARY AND OUTLOOK 



We have studied the relativistic BCS-BEC crossover for zero and nonzero temperatures within a boson-fermion 
model. Variations of this model have previously been used for nonrelativistic systems in order to study cold fermionic 
atoms and high-temperature superconductors. The bosons of the model are bound states of fermion pairs. Conversion 
of two fermions into a boson and vice versa is implemented by requiring chemical equihbrium with respect to this 
process. The crossover is reaHzed by varying an effective coupHng strength x, constructed from the difference between 
the renormalized boson mass mb^r and the boson chemical potential p,b, and the boson-fermion coupling constant g, 
X = —{ml^ — /i^)/(4f/^). In this form, 1/x plays the role of the scattering length, in particular 1/x = ±oo in the 
unitary limit. We have evaluated the model in its simplest form, employing a mean-field approximation. 

An important property of the model is the coexistence of weakly-coupled Cooper pairs with condensed and uncon- 
densed bosonic bound states. In the crossover regime as well as in the BEC regime, strongly-bound molecular Cooper 
pairs exist below and above the critical temperature Tc. Above T^, they are all uncondensed ("preformed" Cooper 
pairs) while below Tc a certain fraction of them forms a Bose-Einstein condensate. In contrast, in the BCS regime, 
pairing and condensation of fermionic degrees of freedom (in the absence of bosons) both set in at Tc. 

Furthermore, we have characterized the onset of nonzero antifermion and antiboson populations during the crossover 
by computing the energy density. The reason for the appearance of antiparticles is the strong decrease of the fermion 
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Figure 10: The phase diagram in the plane of the crossover parameter x and the density difference Sn/n. Shaded areas have 
unstable homogeneous solutions with negative number susceptibility. NS denotes "normal state"; in this region, no solution for 
the gap equation is found. gSF denotes "gapless superfluid"; in this region a stable gapless superfluid state is found with two 
different Fermi surface topologies, divided by the right dashed-dotted line. The labels II, III, and IV refer to the states listed 
in Table U 

chemical potential. While the fermion chemical potential is identical to the Fermi energy in the BCS regime, it reaches 
values well below the fermion mass in the BEC regime. As a consequence, particle and antiparticle excitation energies 
become almost identical and thus antiparticles are present for nonvanishing temperatures. 

Finally, we have extended the model by considering two fermion species with mismatched densities. This case has 
been evaluated for zero temperature. We have found stable gapless superfluids in the crossover region. In contrast to 
nonrelativistic systems, we found no stable homogeneous phase in the far BEC region. Moreover, two different stable 
Fermi surface configurations have been identified. Besides a state with a single effective Fermi surface, also found in 
nonrelativistic systems, we found the possibility of a superfluid phase with Fermi surfaces for particles of species 1 
and anti-particles of species 2. A complete evaluation of the two-species model, including inhomogeneous phases and 
nonzero temperatures, remains to be done in the future. 

The model may be extended in several ways, in order to describe more realistic scenarios, for instance dense quark 
matter in the interior of a compact star. First, one may go beyond the mean fleld approximation, which seems 
particularly interesting in the crossover region, where the validity of this approximation is questionable. Also, one 
may introduce more than two fermion species, accounting for color and flavor degrees of freedom in quark matter. 
Finally, we propose to include chiral condensates into the model. 
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In this appendix, we compute the elements of the number susceptibility matrix as given in Eq. (|55|) . For a compact 
notation we introduce the following abbreviations for integrals containing a (5-function 
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and integrals containing a step function 



e 

0^ - E 



(27r)3 2(6^)3 



(27r)3 4(6^)3 

In this notation the various terms in Eq. (|55|) become 



and 



5A 



dfJ.2/1 



= 02-63. 



dfJ.1/2 

The partial derivative of the gap with respect to the chemical potentials, obtained from the gap equation, is 

9 A _ 2^A/.g2 + 61 T<5f 



1/2 



4(62 - 53) 



Consequently, 



Xii/22 



A^^^ , (2MA/g^+e^^5±)2 

2^2 4(02 - (S3) 



_ A2 (2/iA/g2 + 9i - (?+)(2A2A/.g2 + 61 + 6^) 
X12 - X2i-^ + e2-<53 + 4(62 ~<53) 

In particular, we see that the susceptibility matrix is symmetric. 

We may further evaluate the terms with the (5-functions. For any function /"^(fc) we have 



(A2a) 
(A2b) 

(A3) 
(A4) 

(A5) 

(A6a) 
(A6b) 



E 



rj/ {k)5{5fi-e^) = 



(27r)3 



{e(p+)vpTC+/+(V?T) 



27r2 ^5^2 _ A2 

+ e(p_)v7^C-[e(C-)/+(Vp^) - e(-C-)/-(Vp^)]} , 

with ^± and p± defined in Eq. l(52|) . Consequently, 



<S3 = 



1 Ae((5^t- 


-A) 




- A2 




A) 


27r^ 4(5AiVV 


- A2 


1 A2e(d> 


-A) 


27r^ 4,5^*^(5^^ 


- A2 



[Q{p+)VpTC+ h± + Q{p-)^C~ sgnC-] 

[e(p+)VpTC+ hi + e(p_)VP^C- sgnC-] 
[e(p+)vpTC+ + e(p_)Vprc- sgnC-] , 



where we abbreviated 



h± = Sfi ± ^ySfi^ - A2 . 

For the integrals with the step function we make use of Eqs. ((53|) to find for any function f^ik) 
(Pk 



(A7) 

(A8a) 
(A8b) 
(A8c) 

(A9) 



E/|^me(4-^^)-E/|^/' 



(fc) 



9(<?/^- A) 
27r2 



e(p+) / dkk'f+ik) 



'0 

/p- r^p- 
,2 



QiC-Mp-) dkk'f+ik) + e{-C-Mp-) dkk'rik) 

Jo Jo 



(AlO) 
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We insert Eqs. (|A8P and Eq. I|A10P (the latter with the respective integrand replacing /^(k)) into Eq. (|A6P to evaluate 
the susceptibility matrix x- 
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